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Abstract 

The non-equilibrium dynamics of an amorphous solid is studied with a soft-spin type model. We 
show that the aging behavior in the glassy state follows a modified Kohlrausch- Williams- Watts 
(KWW) form similar to that obtained in Phys. Rev. Lett. 95, 055702 (2005) from analysis of 
the dielectric loss data. The nature of the fluctuation-dissipation theorem (FDT) violation is also 
studied in the time as well as correlation windows. 
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In the glassy state, the liquid behaves like a frozen solid with the motion of its constituent 
particles being localized around randomly distributed sites. Analysis of the dynamics in this 
non equilibrium glassy state reveals a variety of phenomena like aging and memory effects 
Important progress in understanding the non equilibrium dynamics of disordered 
systems has been made in recent years from study of simple mean field spin glass models. In 
the multi-spin interaction models, the non-linearities in the Langevin dynamics give rise^ 
to a ergodic-nonergodic transition. The basic mechanism for this transition is very similar 



to that present in the models for the dynamics of supercooled liquids [4]. The seminal work 
of Ref. [5] dealt with the problem of weak ergodicity breaking^] in a spherical p-spin (p > 2) 
interaction model [7| over the asymptotic time scales. Here the crossing over of the dynamics 
from a regime of time translational invariance to that of aging behavior was demonstrated 
analytically. Low temperature properties of glassy systems, e.g., thermal conductivity and 
specific heat, have also been studied with models 8] for the structural glass in terms of a 
standard Hamiltonian involving spins. In the present paper, we study a soft-spin type model, 
which is defined in terms of the displacements of the particles around a corresponding set 
of random lattice points. We show that the aging behavior in the nonequilibrium glassy 
state follows a modified KWW form similar to that obtained in Ref. |9| from analysis of the 
dielectric loss data for several materials below the glass transition temperature T g . 

We consider a model Hamiltonian, which has a translationally invariant form in terms of 
the displacement variables Ui around an amorphous structure. 

oo 

ff = £E-# ) (*-«') p • (i) 

p=2 v£j 

For the amorphous solid, the interaction matrix J-?' is assumed to be random following a 
gaussian probability distribution of zero mean and variance J£/N. The microscopic basis 
for such a model for an amorphous solid is discussed further below. The time evolution of 
Ui(t) is given by the dissipative Langevin equation, 

r is the bare kinetic coefficient related to the variance of the gaussian white noise £j through 
the fluctuation-dissipation relation < £i{t)£j(tf) > = 2/3~ 1 r 8ij5(t—t / ). z{t) is a Lagrange's 
multiplier used enforcing the constraint N~ l £\ < u?(t) >= 1. This in the present context 
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is equivalent to having a constant Lindemann parameter at a fixed temperature T. The 
simplest form of the nonlinear Langevin equation (j2J) is obtained by keeping in H the con- 
tributions from the p = 2 and 3 terms of the expansion ([I]). In the present work, we focus 
on the corresponding nonlinear model for the dynamics and refer this as the p23 model from 
hereon. 

We formulate a standard Martin-Siggia-Rose (MSR) [k3| field theory to compute the time 
correlation and response functions of the displacement variable Ui. The field iii conjugate to 
Ui is introduced in this regard to average over the gaussian noise £j. The two time correlation 
and response functions are respectively defined as : C(t,t w ) = N' 1 YuiLi ( u i(t) u i{tw)) > and 
R(t,t w ) = N" 1 ^^ (ui(t)ui(t w )) , where the overbars stand for averages over the random 
bonds {J-? '}'s and the angular brackets represent mean over the gaussian white noise £»'s. 
The dynamics of the correlation and response functions are obtained from the equations 



ft rt v 

[d t + z(t)]C(t,t w ) = / dsE(t,s)C(s,t w ) + / dsE'(t,s)R(t w ,s) (3) 

Jo Jo 



[d t + z(t)]R(t,t w ) = 5(t-t w ) + / dsE(t,s)C(s,t w ), (4) 

Jt v 

where we denote E'(t, t') = 25 (t — t') + 5(t, t'). The kernels are obtained from a perturbative 
summation as E(t,t') = z~2p a P CP ~ 1 ( t 7 t '): and ^(t^') = E P (P ~ l)a p C p - 2 (t,t')R(t,t') in 
terms of a set of coupling constants {a p }, which depend on nonlinearities in the dynamic 
equations. For the p23 model, we obtain up to one loop order a 2 = 2((3J 2 ) 2 and = 
ld>((3J?) 2 . The necessary boundary conditions for C and R are respectively chosen as 5] : 
R(t,t~) = 1 and dfCit,^) = ±1. The Lagrange's multiplier z(t), which ensures C(t,t) = 1, 
is obtained as z (t) = 1 + / ds {E(t, s)R(s, t) + £(t, s)C{t, s)}. 

The analysis of the asymptotic dynamics of C(t,t w ) for both t and t w — > oo is divided^] 
into two main regimes. First, for (t — t w )/t — > the time translational invariance (TTI) 
holds. At this stage C and R are related through the fluctuation dissipation theorem (FDT) 
i? : (t) = -Q(t)dtCi(t), where we denote C(t+t w , t w )=Ci(t) and R(t + t w ,t w )=Ri(t). Second, 
for (t — t w )/t ~ 1, ie., for widely separated t and t w there is aging behavior. The correlation 
and response functions, respectively denoted by Ca and R& in this case, are assumed to be 
functions of t w /t = A (0 < A < 1). We define C A (t,t w ) = qC(X) and R A (t,t w ) = t~ l 1l(\). In 
the limit, A — > 1, C(A) — > 1 and 7Z(X) ^ 0. The solutions in the FDT and the aging regimes 



2 



agree if the long time limit of C\(t) is q termed as the non ergodicity parameter (NEP). In 
the FDT regime, both the eqns. ([3]) and (jlj) reduce to a single equation 

(d t + 1) d(t) + [ dsE F (t - s)9 a Ci(s) = Zoo[Ci(t) - 1]. (5) 
Jo 

The kernel H/[Ci] reduce to a^Ci + aaCf in case of the p23 model. Except for the linear term 
on the RHS, eqn. (jSJ) is same as the basic dynamical equation in the self-consistent mode- 
coupling theory of the structural glass. The latter represents the asymptotic dynamics for 
the time correlation of the equilibrium density fluctuations in a supercooled liquid. However, 
in the present case the nontrivial renormalization contribution to the transport coefficient 
comes from the dissipative nonlinearities in (j2J), while in the MCT for compressible liquids 
the relevant nonlinearity is in the reversible pressure term. From the t — > oo limit of eqn. 
(JSJ) we obtain the following relation 



J2a p (p-l)q p - 2 +(l-q) 2 = . (6) 

v 

for the NEP q in terms of the coupling constants a p . 

In the aging regime, the FDT violation is denoted in terms of a parameter m, which is de- 
fined through the relation -Ra(^) = —TnQ(t)d t Cx{T) or equivalently 7£(A) = —mq(d/d\)C(X). 
We obtain analyzing eqns. (j3J) and (|1J) in the aging regime, the following equations for m 
and q, 

At the transition m = 1. The critical coupling constants {a^, otg} for dynamic transition point 
of the p23 model is obtained from the solution of eqns. (jSJ) and ([7j) as = 2/A — 1/Aq 
and ag = 1/Aq, where Xq = 1 — q. The ergodic-nonergodic transition line is given by 



a* 2 = 2^/a^ — 03. This is identical to the line of dynamic transition in the 0i2 mo del [111. Il2 ] 
of the mode coupling theory of structural glass transition. Along the line of transition the 
parameter Ao changes from to 0.5 as the NEP changes from 1 to 0. In the ergodic phase, 
the NEP q = and the FDT holds with m—1. Close to the transition line, the relaxation 
behavior follows several regimes crossing over from power law decay to a final stretched 



exponential form. The corresponding stretching exponent (3^ is approximated 
empirical relation y p a p qp~ 1 ^^ = 1. 
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To get a better understanding of the time scales associated with the aging dynamics in 
the intermediate time regime and the corresponding FDT violation, we solve the eqns. ([3]) 
and flU) numerically. This requires integrating the equations (EJ) and (j3J for both t and £ w 
extending over several time decades. We use the adaptive integration technique 13j], which 
starts with smaller sized grids for integration over shorter time scales of fast relaxation and 
correspondingly increases the step size for longer time scales of slow dynamics. In the ergodic 
state, at long waiting times t w the correlation function approaches its equilibrium value and 
time translational invariance is eventually reached. In the approach to the equilibrium, the 
waiting time dependence of C(t + t w , t w ) is displayed w.r.t. t in fig. CQfor ci2 = 0.82 and a 3 = 
2.02. In the final stage, the decay follows the stretched exponential form exp[— (t/r^ E )^ E ] 
with characteristic relaxation time t^ e and stretching exponent /3^ E . The inset of fig. 
[1] shows f3^ E corresponding to different waiting times t w . At large £ w , it approaches its 
equilibrium value /3 E , which is determined in terms of a 2 and a 3 using the empirical relation 
discussed above. 

In the non-ergodic state, the numerical solution of eqns. (J3J) and (T4J) displays both FDT 
and aging behavior. In figEJ the time dependence of C(t + t w ,t w ) corresponding to a 2 = 0.5 
and a 3 = 6.0 deep in the glassy state are shown for different values of t w . Initially the 
correlation decays from 1 to q and at this stage time translational invariance holds. The 
dynamics is strongly dependent on t w at a later stage. The corresponding correlation and 
response functions for large t w are scaled with the ansatz : C(t + t w , t w ) = C[h(t + t w )/h(t w )] 
where h(t) is a monotonically ascending function of t. The simplest possibility h(t) = P 
is termed as the simple aging and implies C(t,t w ) = C(t/t w ). We adopt here the more 
general formQ, [ijj] h(t) = exp[t 1_K /(l — k)]. The limit k— >0 implies time translational 
invariance while k— >1 represents simple aging. The case < k < 1 is termed as sub-aging. 
The dynamics almost conforms to simple aging behavior as shown in the inset of figj2] in 
which different t w data overlap on a single master curve having k = 0.96. For every t w , 
the correlation C(t + t w ,t w ) decays to zero at sufficiently long t. This is termed as weak 
ergodicity breaking in the aging regime. 

We now focus on the aging time dependence of the relaxation. For a set of £ w 's, the 
fourier transform of the correlation function C(t + t w ,t w ) with respect to t is obtained 
numerically. Since the correlation function in the aging regime is approximately function of 
t/t w (k = .96), the corresponding fourier transform C(u,t w ) is a function of ut w = t w . For 
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comparison with experimental data, we define the response function x w (t w ) = uC(u,t w ). 
The waiting time (f w ) dependences of Xw(^w) f° r different frequencies do not fit with a simple 
stretched exponential form exp[— (t/r) 13 ] with constant r over the whole time range and a 
frequency independent f3. The data is fitted with the modified KWW in a manner similar 
to that of Ref. Q. 

xUQ = [xS - X^ q ] exp [-(t w /r(t w ))^ + X 7 (8) 

where the subscripts "st" and "eq" respectively refer to the limits t w — > and oo for 
The aging time dependence of r is chosen as 

r(t w ) = {r st - r fn } /(f w ) + r fn , (9) 

where r st and Tf n are fit parameters independent of frequency to. The normalized function 
f(s) is chosen to have limiting values 1 and for s — » and oo respectively. In particular, 
we make the choice [id] f(s) = a a /[l + exp {s/r(s)} /3 ] where a = 2@ is a normalization 
constant. Using this form of the r(t w ) we have fitted Xw(tw) for all the different frequencies 
with a single ( frequency independent ) stretching exponent (3. In fig. [HI a scaled plot of the 
different frequency data with respect to t w is displayed. The data sets for all the frequencies 
merge on a single master curve with f3 = 0.55 and is shown as a solid line. For the dielectric 
loss data, Lunkenheimer et. al. in Ref. [9[ use a somewhat different fitting scheme with 
the f(s) in eqn. (jUJ) being a stretched exponential function exp[— (s/r(s))^]. But these 
authors adopt the parametrization of eqn. ([9]) not for the time dependence of the relaxation 
time r(t w ) in the modified KWW formula, but for the corresponding frequency defined as 
u(tyf) = l/{27rr(t w )}. Relaxation data when fitted with this scheme obtain the exponent 
(also frequency independent) f3 = 0.53. In the inset of figlH the r's from both of the above 
described fitting schemes are displayed. The stretching exponent values are close although 
the relaxation time r(t w ) are in fact quite different in the two schemes. 

We now consider the FDT violation in the nonequilibrium state. The FDT is generalized 
in terms of a quantity X(t,f) ( for t > t' ) as k B TR(t,t') = X(t,tf)dC(t,f)/dt'. In the 
limit t,t' — > oo, it is assumed that X(t,t r ) = x[C(t,t')} representing FDT violation in the 
correlation windows rather than time windows. For convenience of discussion, an integrated 
response function F(t,t') is defined 
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F(t,t')= / dsR{t,s) = -— x{C)dC . (10) 
Jt> k B T Jc 

If the FDT holds, x — — 1 and the above relation reduce to ksTF(t) = C(t) — 1, using 
C(t,t) = 1. An effective temperature T e g for the nonequilibrium state is defined in terms of 
the ratio of the fourier transforms, fceTF(u;, t w )/C(u, t w ). If the FDT holds, T e s = 1. Using 
the relation (j7]), we obtain that the choice a<i = 0.5 and 03 = 3.0 in the p23 model makes 



T e fj close to the experimental result of Ref. 17|. More importantly, the time scale of t 



over which the cross over from the FDT to the aging regime occurs according to the present 
model is comparable with experimental observations as shown in fig. HI We also display in 
the inset of this figure, the FDT violation corresponding to the case of fig. [2] as seen from 



the correlation windows. This is similar to results 18] from molecular dynamics simulations 
of the binary Lennard- Jones mixtures. 

The present model for an amorphous solid can be justified from a semi- microscopic basis. 
The potential energy is expressed as a Born von Karman type expansion of the coordinates 
{rj} of the N particles, 

H = ^ JifuiUj + ^ J ijl U i U i U k + •• + G(v,i), (11) 

ij ijk 

where r, = rf + «j. The primes in the summations in the RHS indicate that the terms 
having all the corresponding running indices k etc. being same are absent. In case of 
the amorphous solid, {r 4 } constitute a random structure corresponding to a local minimum 
of potential energy and itj is the displacement of the z-th particle from its parent site. 
The expansion in terms of u^s is valid over the time scale of the structural relaxation. 
The single site potential G(ui) |8| in the RHS of f TTTT) is being included to stabilize the 
system. We will approximate U = J2i<j 4>ij as a sum °f ^ wo body potentials and write 
the potential energy in the translationally invariant form given by eqn. (pQ) by assuming 
Jijl = Sjk + jfk^ki + J/d 8ij an d Jjf = (—l) p jjf etc. For reaching the expression 
(CQ) the coefficients of the single site term G(ui) = J2i{ w 2i u l + w 3i u i + •••} are chosen as 
: W2i = —^2'jJij\ = ~Yl'j k Jijh e ^ c - The semi-microscopic interpretation described 
above is useful in linking the model with thermodynamic parameters [19] . This will test 
further the possibility of using the mode coupling approach to study the complex dynamics 
of the non-equilibrium state of an amorphous solid. CSIR, India is acknowledged for financial 
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FIG. 1: The correlation C(t + t w ,t w ) vs. t for a set of waiting times i w 's in the ergodic phase, 
a>2 = .82 and 03 = 2.0. Inset shows the exponent j3 a for final stretched exponential relaxation w.r.t. 
t w . Dashed line is the corresponding equilibrium value of stretching exponent. 
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FIG. 2: The correlation C(t + t w ,t w ) vs. t for a set of waiting times t w 's in the nonergodic phase, 
a2 = 0.5 and 03 = 6.0. Inset shows scaling of different C(t+t w , i w )'s as a function of h(t+t w )/h(t w ). 
Here h(t) = exp[t 1 ~ K /(l - «)] with « = .96. 
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FIG. 3: Scaled function (xoj(t w ) — xff)/(xS 3 — Xw) vs - *w Different symbols represent data for 
different frequencies u. Solid line is a fit to modified KWW form with r(t w ) and (3 = 0.55. Inset 



shows the variation of the r with i w , following the scheme of (a) Ref. 
(dashed). 
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FIG. 4: Effective temperature T e g (see text) vs. waiting times £ w for u = 10 — 5 (filled circle), 3 x 
10 -4 (open circle), and 5 x 10 -4 (filled square) in units of . Solid line corresponds to best fit curve 
of theoretical results. Experimental data of Ref. [ljj shown with filled triangles. Inset shows 
—F(t + t w ,t w ) vs. C(t + t w ),t w ) corresponding to the data of fig. [2j The slope in the FDT 
violation regime is m = 0.465. 
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